function y=Compute_C_j(a,b)
    s=length(a)-1;
    y=zeros(5,1);
    y(1,1)=sum(a);
    for i=1:4
        for m=0:s
            y(i+1,1)=y(i+1,1)+a(m+1)*(m^i)/factorial(i)-(b(m+1)*(m^(i-1))/factorial(i-1));
        end
    end
end